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The Bethe-Salpeter amplitude is expanded on a hyperspherical basis, thereby reducing the orig- 
inal 4-dimensional integral equation into an infinite set of coupled 1-dimensional ones. It is shown 
that this representation offers a highly accurate method to determine numerically the bound state 
solutions. For generic cases only a few hyperspherical waves are needed to achieve convergence, 
both for the ground state as well as for radially or orbitally excited states. The wave function is 
reconstructed for several cases and in particular it is shown that it becomes independent of the 
■ relative time in the nonrelativistic regime. 

I. INTRODUCTION 

Over the past 40 years the Bethe-Salpeter equation has been the traditional starting point for relativistic field 
theoretical studies of two-body bound states . In view of the numerical complexity, most studies are restricted 
to the ladder approximation. Within this approximation the two-body amplitude obeys a 4-dimensional integral 
t-H \ equation that can easily be transformed into a Fredholm-type eigenvalue problem by a Wick rotation Except 
for the Wick-Cutkosky model |i|,|5| where two scalar particles interact through the exchange of a massless particle, 
very limited analytically obtained information is known about the solutions of the ladder Bethe-Salpeter equation. 

Moreover, despite its fairly transparent structure, it has proven to be a nontrivial numerical problem to find accurate 

OV . . 

solutions of it for a more general type of interaction. 

Many studies have been carried out to determine the Bethe-Salpeter eigenvalues and vertex functions numerically. 
Most efforts start with the reduction of the 4-dimensional integral equation into a 2-dimensional one by introducing 
a partial wave decomposition, thereby exploiting rotational symmetry The 2-dimensional integral equation is 

. then solved by standard discretization methods. Although the resulting matrix equations can be solved numerically, 

•i-H . 

one has often to deal with very large matrices, especially when the spin complication of the external particles has to 
be accounted for. In addition, this happens also when the exchanged mass becomes small. To account properly for 
the long range character of the interaction, a larger number of integration points is needed in that case. Variational 
calculations |^J| have also been done and proven to be a powerful technique. They are, however, harder to apply to 
studies of excited states. 

In this work we consider an alternative method for solving the ladder Bethe-Salpeter equation based on a hyper- 
spherical expansion ]Tc[ | . Instead of doing a standard partial wave decomposition, thereby expanding the amplitude on 
a basis of 0(3) spherical harmonics, we expand the amplitude on an 0(4) basis. As a result the ladder Bethe-Salpeter 
equation is transformed into an infinite set of coupled 1-dimensional integral equations. This representation clearly 
has the advantage that the integral equations are substantially easier and more accurately solvable, provided that 
the number of channels is limited. The idea behind this method is in fact quite old and some approximated, lowest 
order results have already been published a long time ago [ pd| . Here we will systematically improve upon the lowest 
order calculations and solve the full set of integral equations until the desired accuracy has been achieved. For typical 
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calculations only very few hyperspherical waves are needed in practice in order to find a reasonable accuracy. In 
particular when the exchanged mass or the binding energy is sizable as compared to the constituent masses, only 1 
or 2 hyperspherical waves are required to reach convergence. 

In Sect. O we discuss the ladder Bethe-Salpeter equation and introduce the hyperspherical expansion to derive the 



resulting coupled set of 1-dimensional integral equations. In Sect. Ill we present the numerical results. The rate of 



convergence of the eigenvalues and vertex functions is studied and shown to be fast which renders the method an 



efficient one. Finally in Sect. IV we show that the method can also be used to reconstruct the solutions in the case 
of (2+1) dimensions and some concluding remarks are made. Two appendices contain the explicit expressions of the 
kernel of the Bethe-Salpeter equation on the hyperspherical basis, both for (3+1) as well as for (2+1) dimensions. 

II. THE LADDER BETHE-SALPETER EQUATION 

For definiteness, let us consider the bound state problem of two scalar particles with mass m in (d — 1) spatial and 
1 time dimensions. In a relativistic field theory a two-particle bound state with mass M corresponds to a solution of 
the homogeneous Bethe-Salpeter equation (BSE) at total momentum P with P 2 = M 2 . After a Wick rotation [Q, 
the BSE for the corresponding vertex function $(g) has the general structure 

<%) = 7^7 / dV V(q, q')S (q') <&(<?'), (1) 



(2ir) d 

where q is the relative momentum, V (q, q') is the set of all 2-particle irreducible diagrams and S(q) is the free 2-particle 
Green function. In the CM-frame (P — (0, y/s)) we have for S(q) 
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S, ' ?) (q 2 + m 2 - \s) 2 + (q ■ P) 2 ' 

(2) 



(q 2 + m 2 — \s) + q 2 s cos 2 \p q 

with xPq being the angle between the (Euclidean) vectors P and q in d-dimensional momentum space. For a given 
V{q,q'), Eq. ([j]) only supports solutions for values of y/s < 2m that correspond to bound states. 

In this paper we will confine ourselves to the particular example of two scalar particles with mass m, interacting 
through the exchange of another scalar particle with mass /i. In the ladder approximation we have 
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V(q,q')=g 
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(q - q') 2 + / i2 ' 



^ q 2 + q' 2 + ^ 2 -2\q\\q'\c08 lqq r (3) 

with j qq ' the angle between the relative momenta q and q' . Eq. (Q) sums up the set of ladder diagrams depicted in Fig. 
|l|. In this work we focus our attention to the (3+l)-dimensional case. Appendix [b| contains the essential expressions 
for the (2+l)-dimcnsional case and the generalization to other space-time dimensions can be made analogously. 

Usually Eqs. |l]-[|) are solved by noting that for d = 4 they exhibit 0(3)-symmetry and as a result it is profitable 
to expand the amplitudes on a basis of 0(3) spherical harmonics Yi m (9, <fi). This reduces the ladder BSE to a 2- 
dimensional integral equation in the 3-momentum |q| and the relative energy variable q^. Since for s = the BSE 
has an additional 0(4)-symmetry, which as a consequence reduces it to a single 1-dimensional integral equation in 
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that particular case, it is tempting to explore the convergence rate of an expansion in terms of hypersphericals in 4 
dimensions. Such a method should in particular work very well for the case of strong binding. Starting from Eq. (|l]) 
we may expand the vertex function in appropriate 4-dimensional spherical harmonics. From Eq. (§ it is clear 
that for s ^ the BSE-kernel is not 0(4) invariant. Hence the expansion in hyperspherical harmonics will not lead 
to a complete separation of the generalized angular variables. 

In (3+1) dimensions the suitable set of spherical harmonics is given by: 

Z klm {x,e^)=X kl ( X )Y lm {9^) 1 (4) 

with 



Xu{x) - J— [ (fc ^ + 1) J , W X CJ+Kcobx), (5) 



iwm) - (-D-y^fr^l pr{cos9) eim *- (6) 

Here C l k t_\(x) are the Gegenbauer polynomials of degree k — I and order I + 1 and P{ n {x) are the associated Legendre 
functions of the first kind of degree I and order m For a fixed value of k, the index / assumes integer values 

from to k, while m goes from — I to I. Standard properties of these functions allow one to prove completeness and 
orthogonality for the set Zki m (Xi ®-> ( t ) ) with respect to the scalar product at hand 

d(f> dOsxnO dxsin 2 xZklm(x,0, 4>) z t'i'm'(x, 0, 4>) = Skk'5u>5 mm ,. (7) 
o Jo Jo 

By expanding <£>(g) and V(q,q') in Eq. (Q) as 

<%) = ^2$kl m (\q\)Zklm(X,0,<P) (8) 

klm 

V(q,q') = Y, V k(\ll\l'\) C l( C0S ^') ( 9 ) 

k 

and using the addition theorem for the 4-dimensional hyperspherical harmonics 

c£(cos 7 ) = ^£ Wx,M)^(x',0', A (10) 

lm 

where 7 is the angle between the two unit vectors (x, 9, tp) and (%', 9' ', </>'), we obtain after some straightforward algebra 
an infinite set of coupled 1-dimensional integral equations: 

<w«) - 8 ^ 2(fc + 1) E y o d <?v 3 ^fe?') 4 fe '(9') ^w^O- (11) 

Here the matrix S l kk , (q) is given by 



SjU'(<?)= / dxsin 2 x^(x)^i(x)^(l<7l,cosx). (12) 

r our particular case the mat 





For our particular case the matrix-elements 14 (q, q') and S l kk , (q) can explicitly be computed and are given in Appendix 



Eq. (11) exhibits several symmetries that allow for a physical interpretation and which lead to simplifications in the 



numerical treatment. Since Vk(q, q') and S l kk ,(q) are independent of the magnetic quantum number m, the solutions of 
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( |Tl| ) do not depend on it. Hence there is a (21 + l)-fold degeneracy present which is a manifestation of the cylindrical 
symmetry. Furthermore, Eq. (|ll]) is diagonal in the orbital quantum number I which implies the conservation of 
spatial angular momentum, in contrast to the nonconservation of the generalized angular momentum connected to the 
quantum number k. Considering the behavior under parity transformations x — > —x of the Gegenbauer polynomials: 
Cj^](— x) = (— l) k ~ l C k ^} l (x) leads to the conclusion that S l kk ,(q) is only nonzero if k + k' is even. As a result only 
states of equal parity are coupled to each other in ( |ll| ) and consequently parity is conserved. Finally, one can show 
that S l kk ,(q) vanishes unless k > I and k' > I. Hence, for a fixed value of I, the quantum numbers of the generalized 
angular momentum assume values k — I + 2 j and k' = I + 2j', with j 6 0, 1, 2, . . .. 



III. NUMERICAL RESULTS 



Having transformed the original 4-dimensional integral equation (|l|) into an infinite set of coupled 1-dimensional 
ones, our aim is now to find those values of s < Am 2 for which Eq. ([ll]) allows solutions. As noted these would 
correspond to bound states with mass M = y/s. In particular it is interesting to investigate how many hyperspherical 
waves are required in order to achieve a certain accuracy and to see how this number depends on the properties of 
the particular state one is considering. 
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A. Numerical Implementation 

In order to evaluate the integral in Eq. (pj]), we project the infinite interval q = [0, oo) on to the finite unit interval 
z = [0, 1] through the transformation: q —> z/(l — z). The integral can then simply be computed by Gauss-Legendre 
quadratures Pq| . 

N 

dxf(x)= lim y~]wif(xi), (13) 

N—>oo * — ' 
i=l 

with weights Wi and abscissas Xi. Next we truncate the infinite sum over k' in Eq. (|ll]) at a finite value of k' = fc max . 
As a result we have replaced the kernel of the integral equation by an approximated finite matrix K of dimension 
(N ■ fc max ) x (N ■ fc max ) of which we seek the eigenvalues. After having computed det(l — K) numerically, we apply 
Newton's method |l3| ] in order to find the roots of det(l — K) as a function of y/s for a fixed value of g or vice versa. 
Typical values of N = 60 have been used in the actual calculations, except for the case of smallest fi, where N = 200 
was needed to get the same relative accuracy as the other ones. Finally we checked for a large number of states that 
we obtained identical results as zur Linden fll4|. 



B. Convergence of the Eigenvalues 



Since the BSE is only 0(4)-invariant for y/s = 0, the rate of convergence of the hyperspherical method is expected 
to show a significant dependence on the value of the bound state mass y/s = M. We study the convergence rate by 
keeping M fixed at different values and determining the value of the dimensionless coupling constant A = g 2 /Aitm 2 
that yields a bound state of these masses. 
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In Table | we present a comparison of calculations for different values of M and /x. For the calculations we consider 
the ground state where the radial quantum number n — 1 and the system is in an s-state (I — 0). We indeed find 
for a given /i that a larger value of A; max is required for increasing M in order to obtain a certain accuracy. Except 
for the smallest value of /Lt/m = 0.001 and the weak binding case of M/m = 1.999, it is remarkable to see how fast 
the method converges. The latter exceptional case corresponds to the nonrelativistic situation where the binding is 
small as compared to the rest masses of the constituents. A reasonable relative accuracy of 10 -4 can rather easily be 
obtained by considering only a few hyperspherical waves. The convergence is extremely good for moderate values of 
fj,/m ~ 1 where only two hyperspherical waves suffice for the desired accuracy. Comparing the predictions for the case 
of /i/m = 0.001 with those of the Wick-Cutkosky model we find that they agree within 4 digits, except for the weak 
binding situation M/m — 1.999. The poorer correspondence for the latter case is understandable since the exchange 
mass, which vanishes in the Wick-Cutkosky case, is of the order of the binding energy of the system here. Hence its 
influence on the long range structure of the bound state is expected to be significant, while it is negligible for the 
other situations. 

Next we wish to study how the above pattern depends on the orbital quantum number I. For this investigation 
we keep fi fixed at fi/m = 0.1 and we vary I at different values of the bound state mass M. The results of this 
calculation is presented in Table ||. From this table we see, for each value of M, that the rate of convergence is 
virtually independent of the orbital quantum number I. 



Finally, solutions of radially excited states have been studied. Some results are shown in Table [II. Also in this 
case we find that the hyperspherical method works very well. The convergence rate is found to be only very weakly 
dependent on the radial quantum number n. 

From the above investigations, we may conclude that for the computation of the spectrum, the rate of convergence 
of the hyperspherical method is, as expected, dependent on the bound state mass M and on the value of the exchanged 
mass fi. The convergence rate is in general very fast except for small values of the exchanged mass \i. This is not 
surprising since again we expect that the long range character of the force for fj, — ► should be reflected in the rate 
of convergence. In this connection, it should be noted that keeping a finite number of the hyperspherical waves, the 
resulting kernel of the BSE is compact even in the limit fi — > 0. For performing the integrals in an accurate way, more 
integration points are needed when fi becomes small. This may be attributed to the appearance of a cusp in Vk(q, q') 
at q = q' when fi is identically 0. 



C. Hyperspherica[ Components of the Wave Functions 

Clearly, the BSE (jl|) with kernel K has formally the structure of a generalized eigenvalue problem 

A$ = K$. (14) 

After discretization of the integral in the BSE we may find for a given energy ^fs the eigenvalues and eigenvectors of 
the matrix K using standard procedures. The physical solutions can simply be obtained by varying s such that we 
get an eigenvalue A = 1. Once the eigenvectors have been determined of Eq. (|lj) we immediately obtain the wave 
function ^>(q) through the relation 

(?) = (15) 
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In the actual reconstruction of the ground state vertex function we have used the so-called power method p5[ , 
which is considerably faster than the procedures for obtaining general eigenvectors. It simply consists of the iteration 
procedure 

$(«+i) =K$ (n)_ ( 16 ) 

One can readily show, that the limiting function $( Tl ^°°) becomes proportional to the eigenvector with largest eigen- 
value A. This eigenvalue is simply given by A = lirrin^oo /$("). To accelerate the iterative procedure we have 
chosen a starting vector <I> < ' - ) which incorporates the asymptotic behavior of the different hyperspherical waves 

**)(«) 9 ~° <A (17) 



^°V*- a , (18) 



More specifically, we choose $^ (g) according to: 



•Sf(.) = iT7 $^- ("» 

The method was used as an additional check on all numerically obtained eigenvalues. 

Obviously the eigenvectors can be obtained in this way up to a normalization. We have normalized the wave 
functions according to the well known BSE normalization condition given by H 

(2tt)4 (2tt)4 W dP M 



[(27r) 4 <5 4 (g - g')^) + Vfa q')] **(q')=2P, (20) 



Inserting the expansion (^) in (^TJ) we find in the CM- frame for each partial wave /: 

1 f°° 

2(2^4 E/ Q dS^^wC?) - <W (g 2 + m 2 - is)] *fc/,(g)=l (21) 



A' A' 

where we have defined: 



* kl (q) = Y,S l kk ,(q)®kn(q) (22) 



^L' = / d X sin 2 X X kl ( X ) X kn ( X ) cos 2 X 
Jo 



(23) 



In Appendix ^| the explicit expression for A l kk , is given. 

In Fig. |^ we show as typical examples the rate of convergence of the hyperspherical method for the vertex function 
<3> ( | q| , (74) of the ground state for two values of the relative energy in the case of a bound state mass M = 1.9m. 
For the whole |q| — 94 region we need only 4 hyperspherical waves for this case to get a reasonably converged result 
for the vertex function. As we can see from Fig. ||, the |q|- and g4-dependence of the calculated <E> are smooth and 
the range of fall-off in these two variables are similar. 

The rapid convergence of the hyperspherical method is related to the magnitude of the hyperspherical components 
in the wave function. In Fig. ^|we display for I — and fi/m = 0.1 the components of the ground state wave function 
for the cases M/m = 1 and 1.9. The corresponding coupling constants are given by A = 20.678 respectively A = 5.227. 
For the latter case we see that the contributions of higher (k > 2) hyperspherical waves are indeed more significant 
than for the former situation, in agreement with the slower rate of convergence of the eigenvalue. Moreover, from the 
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figure we see that the range of the fall-off in q is very different for the two cases. This is predominantly caused by the s- 
dependence of the free 2-particle Green function and related to the size of the composite system. The hyperspherical 
components of the corresponding vertex functions exhibit a more comparable range and the non-leading ones are 
relatively more suppressed in magnitude. 

In the weak binding regime we expect both particles to be nearly on their mass-shell, which results for the wave 
function ^ in a strongly peaked structure in the q4-variable as compared to the typical fall-off in the |q|-variable. As a 
consequence we expect that the wave function in the configuration space will exhibit an independence on the relative 
time t in the nonrelativistic limit. We have studied this transition region. The wave functions in configuration space 
can be obtained from the above ones via Fourier transformation. With the help of the representation of the complex 
exponential 

Jk+i(qr) 

klm 

we find that for a fixed /, the configuration space wave functions are given by 

*,(M,t) = ^fE^^T^) / <M +1 ( q V^T¥) *«(*) (25) 

The integral Eq. (Pq) was performed by Gauss-Legendre integration. In order to compute the wave functions 1 QE r / ( | r | , 0) 
at relative time zero, one needs the value of X k i(^ir) 



= 4tt 2 V Jk+m ' 1 Z klm ( Xq , 6 q , <t> q )Z*( Xr , r , &•), (24) 



cr-;(o)=(-i)»*-"(j5:;;j - I «'>)=H)'« 1 ( T ^wif (26) 

whereas the determination of the wave function ^/(O, t) at zero relative distance requires X k i(G) 

C l k t)(l)=Si (k+l) =► X kl (0)=5 l0 ^(k+1), (27) 

We display in Fig. ^| the computed configuration space components 

1 f 00 

*w(r) = —t- / dq q 2 J k+1 {qr) * kl {q) (28) 



47r 2 r 

of the ground state, again for M/m =1.0 and 1.9 with /i/m = 0.1. Comparing this with Fig. ||| we see that the 
pronounced structure of the higher components is relatively more suppressed and its strength is spread out to larger 
values of r. The asymmetric behaviour of the configuration space wave function ^/ from Eq. |2^ as a function of |r| 
and t can clearly be seen in Fig. || for the case of M/m— 1.9, 1=0. Its fall-off in t is twice as slow as compared to that 
in |r|. 

Finally, it is of some interest to study the role of the relative time in the transition region to weak binding. In Fig. 
[?] we show three graphs with M/m = 1, 1.9 and 1.999 where we compare the dependence on the Euclidean relative 
time t to that of the relative distance |r|. We expect that in going to the region where the binding energy becomes 
negligible as compared to the total mass of the system, to recover the Schrodinger wave functions and to lose the 
dependence on the relative time which is inherently absent in a nonrelativistic description. For M/m — ► 2 we indeed 
observe that the dependence of \&o(l r l; t) on t becomes substantially weaker as compared to that on r, i.e., the wave 
functions clearly show the expected independence on the relative time t. 
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IV. CONCLUDING REMARKS 



We have investigated the hypersphcrical expansion for solving the ladder Bethe-Salpeter equation in (3+1) dimen- 
sions. It was found that it constitutes a powerful and efficient method for obtaining an accurate determination of 
the eigenvalues as well as the eigenfunctions. This is true both for the ground state as well as for excited states. We 
stress that all calculations presented here took at most one minute of CPU-time on a HP 735- workstation. Except 
for the limiting cases of vanishing exchange mass or small binding energy, only a few hyperspherical waves suffice for 
computing the solutions to the ladder Bethe-Salpeter equation within an accuracy that is sufficient for most practical 
purposes. Furthermore, it was found that the eigenfunctions could be reconstructed in a straightforward fashion. 

The method can equally well be generalized to the case of d-dimension. In particular, we may consider the case of 
2 spatial and 1 time dimensions. The explicit form of the kernel of the BSE for this case is given in Appendix |b| In 
Table |fv| we show the converged results for the dimensionless coupling constant A = g 2 /Anm for fixed values of the 
bound state mass and \x — ► 0. Similar rates of convergence properties are found as for the (3+1) dimensional case, i.e., 
for weak binding and smaller exchange mass \x more terms in the hyperspherical expansion are needed. From Table ^ 
we see that the limit of fj, — > is smooth and no logarithmic confinement is found, as would have been expected from 
a nonrelativistic analysis. Such a logarithmic behavior can readily be seen using the frequently adopted procedure to 
obtain a static nonrelativistic potential Unr(|i"|) from a relativistic interaction: 



VMM)= / dtV(\r\,t). (29) 
For the Yukawa interaction in (2+1) dimensions we find consequently: 

V N R,Yukawa(|r|) = -f-tf (/i|r|). (30) 

The asymptotic behavior of the Bessel function for small arguments Kq(z) ~ — log(|z), suggests clearly that the 



nonrelativistic Coulomb interaction in (2+1) dimensions is logarithmically confining. According to table IV we find 
however the limiting behavior of fi — > to be smooth for the relativistic ladder BSE. Hence no confinement is seen in 
that limit. This paradoxical situation was investigated in [|l6| where it is found that a nonuniformity exists between 
the nonrelativistic limit to — > oo and the Coulomb limit fj, — ► 0. 

Given the very encouraging results we found in this study of scalar particles, we expect the hyperspherical expansion 
to be a powerful and accurate method to solve more realistic physical models such the ones that include spin degrees 
of freedom. In particular, we do expect that it can equally well be applied with good convergence properties for cases 
like the spinor-spinor interaction. Moreover, it should be interesting to investigate the applicability of the method to 
the scattering case. 

APPENDIX A: HYPERSPHERICAL PROJECTION OF THE BSE 3+ i 

In this appendix we give the explicit formulae for the kernel of the (3+1) dimensional BSE that occur in the 
hyperspherical expansion. Most expressions are obtained with the help of standard properties of special functions, as 



can be found, e.g., in Ref. [|12|,[17 18 



The projection of the interaction on the hyperspherical basis can be performed analytically 
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2 

V k (q,q') = g 2 - / dxsin 2 xCl(cosx)V{q,q',cosx) 
~ Jo 



(A+ + A_r 



A+- A_ 
A+ + A_ 



where A± = y/(q± q') 2 + /j, 2 . 

Substituting Eq. (||) in Eq. ( [l2| ) and introducing x — cos ~fp q , we may rewrite S l kk ,{q) as 



sL'(q) 



2 l+1 u 2 (k + i){k-i)\ (k> + i)(k> - iy. r 1 



After some tedious algebra this reduces to 



dx (1 - x 2 ) l+1/2 CftWCfoWSiq, x). 



, ,x_2>/_U! / (fc + l)(fc-Q! [(k> + l)(k'-l)\ - {¥ -\)«Hl I z^+ic^dz^ (i- 
S kk ,(q)-—^^ _____ y ______ e (! + _)* (i_)Q fc _ +I+ ,(u 



with 



g 2 + m 2 - i-s 



(Al) 



(A2) 



(A3) 



(A4) 



and fc + = max(fc, fc') and k_ = mm(k, k'). Despite the appearance of the imaginary unit 'i' in Eq. (|A3|), S l kk ,(q) is 



always real. For I — a more transparent form can be found from (A3). 



s° kk> (q)-(-ir + - k - )/2 



sq 2 z\/l + z 2 



(v/Th - -zy ++k +2 + (-i) fc - (yr+^+z} 



feH fe_ 



(A5) 



The matrix-element A l kk , , defined in Eq. ( p3[ ) as the expectation value of the operator cos 2 x between two spherical 
harmonics, can be computed with the help of the recurrence relation for Gegenbauer polynomials 



2(n + m + l)xC™ +1 (x) = (n + 2m)C™{x) + (n + 2)C™ +2 (x), 



which yields the following expression for A l kk , 



A 1 - i 

A kk' — i 



<>k,k' 



(k-l)(k + l + l) (k-l + l)(k + l + 2) 



k(k + l) 



+Sk.k' 



k,k'+2- 



(fc + l)(fc + 2) 

{k-l){k + l + l)(k' - I + l)(k' + 1 + 2) 



fc(fc + l)(fc' + l)(fc' + 2) 



k' 



(A6) 



(A7) 



APPENDIX B: HYPERSPHERICAL PROJECTION OF THE BSE 2+ i 

In (2+1) dimensions the spherical harmonics are the Yi m (8, </>) functions. The projection of the ladder BSE in this 
case closely follows that of the (3+1) dimensional case. After expanding the vertex function and the interaction on 
the hyperspherical basis 

$(?) = _*Jm(?)lSm(M), (Bl) 
lm 

V{q,q') = Y, V MM'\)Pl{^l^l (B2) 
i 

and applying the appropriate addition theorem 
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A (cos 7) = 2^iE^(^^'^') 



(B3) 



we find 



1 f°° 

In Eq. (B4) the matrix Sfi)(q) is defined analogously to Eq. ( |l2| ) 

SJP(?)= / d6»sin6'P; m (cos6')P^(cos6l)S*(|(?|,cos6l). 
Jo 



(B4) 



(B5) 



Similar considerations concerning conservation of angular momentum (with quantum number to) and parity as have 
been made earlier for the (3+1) dimensional case, apply here as well. 
For the one meson exchange diagram we have 

21 + 1 



Vi(q,q') = 



d6 sin 6»P/(cos 0)V(q, q' , cos ( 



2 2/ + 1 „ 

9 -^rrQi 



2qq> 



q 2 + q' 2 + M 2 
2^ 



Moreover, the free 2-particle Green function with to = 0, S^, (q), can be evaluated analytically: 

yw+iywTT i 



47ri 



sq^z 



(B6) 



(B7) 



where z is defined in Eq. (A4), while ij- are defined according to l + = max(Z,Z') and Z_ = min(Z,Z'). Unfortunately, 
we were unable to find an analytical expression for S^,(q) for arbitrary to. For to ^ 0, the integral in Eq. (B5) was 
performed numerically. 
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TABLE I. Convergence of hyperspherical eigenvalues for A = g 2 /4ivrn 2 as a function of pt/m for n = 1 and I = 0. The 
numerical errors are at most one unit in the last decimal. When no further values are given for increasing fc max , the final 
value was found not to change anymore within the desired accuracy when fc max was increased. The case of fi/m — is the 
Wick-Cutkosky prediction. 



H/m 


^max 


M/m = 


M/m — 1 


M/m — 1.9 


M/m = 1.999 







25.13 


20.01 


4.483 


0.2854 


0.001 





25.13 


20.02 


4.675 


0.3484 




2 




20.01 


4.503 


0.3099 




4 






4.486 


0.2992 




6 

8 






4.483 


0.2947 
0.2924 




10 








0.2912 




12 








0.2904 




14 








0.2899 




16 








0.2896 




18 








0.2894 




20 








0.2893 




22 








0.2892 




24 








0.2891 




26 








0.2890 


0.100 





25.80 


20.69 


5.374 


1.0733 




2 




20.68 


5.237 


1.0463 




4 






5.228 


1.0442 




6 

8 






5.227 


1.0439 
1.0439 




10 








1.0438 


1.000 



2 


42.96 


36.94 


17.28 
17.23 


10.272 
10.251 
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TABLE II. Convergence of hyperspherical eigenvalues for A = g 2 /Aitm 2 as a function of the orbital quantum number I for 
n = 1 and fi/m =0.1. The numerical errors are at most one unit in the last decimal. When no further values are given for 
increasing fcmax, the final value was found not to change anymore within the desired accuracy when fcmax was increased. 



1 


^max 


M/m = 


M/m = 1 


M/m = 1.9 


M/m = 1.999 








25.802 


20.691 


5.3743 


1.0733 




2 




20.678 


5.2369 


1.0463 




4 






5.2279 


1.0442 




6 






5.2271 


1.0439 




8 






5.2270 


1.0439 




10 








1.0438 


1 


1 


160.44 


127.59 


33.284 


10.086 




3 




127.48 


32.094 


9.5620 




5 






32.000 


9.5043 




7 






31.991 


9.4957 




9 






31.990 


9.4942 




11 








9.4940 




13 








9.4939 


2 


2 


668.88 


531.46 


143.58 


55.385 




4 




531.01 


137.99 


51.611 




6 






137.51 


51.123 




8 






137.46 


51.042 




10 






137.45 


51.027 




12 








51.024 




14 








51.023 


3 


3 


2332.0 


1857.7 


525.81 


237.45 




5 




1856.2 


505.29 


219.63 




7 






503.44 


217.12 




9 






503.25 


216.67 




11 






503.23 


216.59 




13 








216.57 
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TABLE III. Convergence of hyperspherical eigenvalues for A = g 2 /4imi 2 as a function of the radial excitation quantum 
number n for I — and \xjrn = 0.1. The numerical errors are at most one unit in the last decimal. When no further values 
are given for increasing fc m „, the final value was found not to change anymore within the desired accuracy when fc max was 
increased. 
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64.387 


16.387 


4.2578 




6 






16.384 


4.2521 




8 






16.383 


4.2510 




10 








4.2508 


3 





169.55 


135.84 


38.053 


12.111 




2 




134.77 


35.557 


10.874 




4 




134.74 


35.257 


10.715 




6 






35.221 


10.690 




8 






35.217 


10.686 




10 






35.216 


10.685 
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TABLE IV. Behavior of the eigenvalues of the BSE in (2+1) dimensions as /i — > 0. The calculations are for the ground state. 



ix, jm. 


M/m = 1.0 


M/m = 1.9 


n i nnn 




U.ooy / 


0.0300 


2.477 


0.2232 


0.0100 


2.366 


0.1902 


0.0030 


2.328 


0.1785 


0.0010 


2.314 


0.1750 


0.0003 


2.308 


0.1738 


0.0001 


2.305 


0.1732 
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FIG. 1. General structure of the ladder diagrams. 
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FIG. 5. Hyperspherical components of the configuration space wave function for fj,/m = 0.1 and I = 0. The upper graph 
shows the components ^ko{r) for M/m = 1 while the lower one contains the solutions for M/m = 1.9. 
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FIG. 7. Comparison of the relative time zero wave function \I>o(|r|,£ = 0) with the relative distance zero wave function 
*()(r = 0, t) for fj,/m = 0.1 and 1 = 0. In the nonrelativistic limit M/m — > 2, *o(|r|, t) is seen to become less dependent on t. 
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